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Abstract 

We report the study of graphene devices in Hall-bar geometry, gated with a polymer electrolyte. 
High densities of 6 xlO^^/cm^ are consistently reached, significantly higher than with conventional 
back-gating. The mobility follows an inverse dependence on density, which can be correlated to a 
dominant scattering from weak scatterers. Furthermore, our measurements show a Bloch-Griineisen 
regime until 100 K (at 6.2 xlO^'^/cm^), consistent with an increase of the density. Ubiquitous in 
our experiments is a small upturn in resistivity around 3 xlO^^/cm^, whose origin is discussed. 
We identify two potential causes for the upturn: the renormalization of Fermi velocity and an 
electrochemically-enhanced scattering rate. 
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I. INTRODUCTION 



Since its first exfoliation from graphite in 2004 [T], graphene transport properties have 
mainly been studied in the vicinity of the Dirac point, where the dispersion relation is linear 
and the electrons behave as massless Dirac particles P-f4J. For technical reasons, the elec- 
trical properties of graphene have rarely been measured at densities beyond 10^^/cm^. But 
the physics of graphene may well be as exciting at high charge carrier densities as it is in the 
vicinity of the Dirac point. As the chemical potential is shifted away from the Dirac point, 
the description of electrons as massless Dirac particles becomes less valid and corrections 
are needed to describe the physics E]. Besides, recent angle resolved photoemission spec- 
troscopy (ARPES) experiments [5] show that potassium- and calcium-doped graphene have 
extended van Hove singularities (VHS), a feature also present in the cuprate energy-bands 
and suspected by some to be responsible for their high-Tc superconducting transitions [5l [7]. 
However, the extended VHS of cuprates are easily accessible, whereas their graphene coun- 
terparts lie at ~ 2-3 eV above the Dirac points, corresponding to electron densities greater 
than 2 x 10^ Vcm^. 

Experimentally, the realization of high carrier densities in graphene devices is limited by 
the requirement of thin dielectrics with high capacitance. These materials are however prone 
to dielectric breakdown at gate voltages required for achieving high doping. In addition the 
growth and identification of graphene on various substrates remains challenging [S]. The 
conventional Si02 back-gate, while being suitable for identifying graphene flakes, cannot 
lead to carrier densities greater than 10^'^/cm^ in graphene. The use of high-K dielectrics 
has also been considered for achieving high doping [9], although this approach has been less 
successful. The present work uses a polymer electrolyte gate to achieve high-doping. When 
a potential difference is applied between two electrodes in an electrochemical cell, the ions 
move in the polymer matrix according to their charge polarity and accumulate to form an 
electric-double layer at the electrode interface. Such nanometer-size gate has a very high 
capacitance and can induce counter charges of equivalent density on graphene. Polymer 
electrolyte top-gating has been previously used to demonstrate the sensitivity of the Raman 
spectrum to high carrier densities in graphene [TOl [TT] . 
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FIG. 1: Optical image and schematic of graphene device in Hall-bar configuration, coated with 
polymer electrolyte [S = Source, D = Drain, G = Polymer Electrolyte Gate]. Scale bar: 10 fim. 

II. POLYMER ELECTROLYTE GATING 

In this work, we study the electronic properties of graphene Hall devices gated with a poly- 
mer electrolyte and track the deviations from Dirac physics through density and temperature 
dependent transport measurements. The Hall measurements demonstrate the effectiveness 
of the electrolyte system in realizing high carrier densities in graphene. From transport mea- 
surements, we evaluate the relative contributions to graphene resistivity induced by different 
scattering mechanisms. Monolayer graphene flakes are prepared by mechanical exfoliation 
on Si/Si02 substrates. The device measurements are performed on standard Hall bar and 
four-terminal structures fabricated with electron-beam lithography. A schematic of the de- 
vice is shown in fig. 1. In addition to the electrodes on the graphene flake, Au/Cr electrodes 
of large surface area are also patterned within the plane of the device structure at few mi- 
cron separations. While several designs for polymer electrolyte gating rely on evaporation 
of top gate-electrode or insertion of a Ft or Au wire in the polymer matrix [12], the 
present design does not require positioning of the top-contact. The in-plane gate electrode 
can be simultaneously patterned lithographically along with the graphene contacts. Migra- 
tion of metal atoms into the polymer matrix may happen in case of evaporated top contact 
and this contamination is also prevented. The polymer electrolyte, an aqueous dispersion 
of polyethylene oxide (FEO) and lithium perchlorate is then drop cast on the device and 
bake-dried. 

The graphene resistance is measured at low frequencies (13 Hz) in the four-terminal 
configuration under low vacuum conditions. A plot showing the modulation of graphene 
resistivity (p) with applied polymer electrolyte gate voltage (Vg) at room temperature is 
shown in fig. 2(a). Due to the large interfacial capacitance arising from a nearby layer of 



counter ions, it is possible to obtain a large and reversible modulation in graphene resistance 
with the application of small voltages. The measurements are restricted to a maximum gate 
leakage of ~ 1 nA. At high gate voltages (or gate leakage currents) the devices show a 
breakdown due to electrochemical reactions. Since the polymer is hygroscopic, the presence 
of adsorbed residual water contributes importantly to this leak [13j, but at the same time 
allows a better ionic mobility. Significantly lower leakage current is observed, when the 
device is cooled below the ice-point of water. In addition, we note that the sweep rate of 
gate voltage must be slow enough to allow equilibration of the ion double layer atop graphene 
and measure a stabilized value of resistivity. The typical mobility of our pristine graphene 
samples at low doping is in the range 4000-7000 cm^/Vs. Upon addition of the polymer 
electrolyte, the mobility of graphene remains larger than 3000 cm^/Vs at ri ~ 10^^/cm^. The 
slope of graphene resistance dR/dVg (measured at half the value of maximum resistance) 
gated with silicon-oxide back-gate is typically 150 This slope is enhanced significantly 

to 3500 Qy~^ when the polymer electrolyte gate is used. Upon sweeping the electrolyte 
gate voltage, a typical on-off ratio of 30-40 is obtained. The sharp resistance slope and high 
on-off ratio value are indicative of high-doping in graphene. At zero gate-voltage, graphene 
is found to be in a highly electron-doped low-resistance state and the charge neutrality point 
is shifted by -3 to -5 V. Such doping may be attributed to a higher concentration of Li+ ions 
adsorbed in the vicinity of graphene, since the graphene has small hole-doping prior to the 
coating of polymer electrolyte [see fig 2(a)]. The G-band Raman peak for graphene shows 
a shift of 6-7 cirT^ upon addition of the polymer electrolyte as well as a reduction in full- 
width at half-maximum (FWHM) [see fig. 2(b)], which further supports the electron-doping 
of graphene [TO] . 

Ubiquitous in our measurements is a small upturn in resistivity observed at high gate 
voltages [see inset of Fig. 2(a)]. This upturn is consistently observed in 6 graphene devices 
on 5 different wafers and across several sweeps for the same sample. To characterize the 
nature of transport at high-doping and examine the contributions to graphene resistivity, 
we performed Hall measurements. The room temperature resistivity p and Hall mobility 
Ph of graphene devices are plotted as a function of the carrier density in fig. 3(a) and 3(b) 
respectively, for two devices. The mobility shows a continuous decrease between 1 x 10^^/ cm? 
and 6 xlO^^/cm^ , /i ~ indicating that p approaches a saturation value. As in the 
low density regime {n < 10^'^/cm^ ), the factors determining the total resistivity include 



4 



(a) 



a 

2 

DC 









1 ■ 1 ' 1 ■ 1 


- 




210 








1 ^ 

1 a: 


200 










190 












12 3 4 6 6 








gate * ' 








1 ' 

~ , . J. . 




1 . 1 







(b) 



-6 -4 -2 2 4 

V , (V) 

gate ^ ' 



3 

re 
> 

c 



c 

to 

E 

(0 



— I — I — I — I 1 — I 1 — I — i— 

» polymer-electrolyte coated graphene 

4 pristine graphene * 

9 3 



_i I I I I L 



1560 1570 1580 1590 1600 1610 
Raman Shift (cm' ) 

FIG. 2: (a) Resistance vs. polymer gate voltage for sample 1 [Inset: R vs. Vgate in the low 
resistance region, showing an upturn in the device resistance] (b) G-band Raman-shift for pristine 
graphene (red) and polymer-electrolyte coated graphene (black). 



charged impurities (both from underlying substrate and from electrolyte ions), defects on 
the graphene lattice and phonons [U [HHTTj. While these contributions have been examined 
at low densities, their relative contributions at high densities can be significantly different 
due to high screening from carriers in graphene. 



III. CONTRIBUTIONS TO GRAPHENE RESISTIVITY 

In this experiment, graphene is sandwiched between the Si02 substrate below and the 
polymer layer above. The number of electrolyte ions in the vicinity of the graphene 
sheet increases with the electron density, unlike the number of charged impurities at the 
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5*^02 /graphene interface which remains constant. To estimate the contribution of the 
ions to the total resistivity of graphene, it is necessary to know their distribution in the 
vicinity of the graphene sheet, which is hard to obtain experimentally. Theoretically, the 
Poisson-Boltzmann equation is often used to describe the ion distribution in electrolyte 
systems p^. However, the concentration of ions estimated from this model diverges at the 
graphene /polymer interface while the concentration of ions is limited by the finite ionic 
radius, the space occupied by the polymer, and the formation of electrolyte-polymer com- 
plex. To take this into account, modified Poisson-Boltzmann equations are generally applied 
and/or cutoff concentrations Cmax introduced [ISj. Following the latter approach, we mod- 
eled the ion-induced resistivity of graphene [see Appendix] . With a polymer packing density 
/ < 80% and an electrolyte ion effective radius around 1 nm, Cmax takes values between 
and 5xl0^^/m^. The polymer dielectric constant is e ~ 5 [lOJ. The concentration 
of ions in the bulk polymer matrix is estimated to be about 5xl0^^/m^. The gate voltage 
dependence of carrier density is plotted in Fig. 3(d). This can be used to experimentally 
estimate the total gate capacitance (polymer capacitance and quantum capacitance in se- 
ries), which is of the order of ~ l^F/civ? [see fig. 3(d)]. Second, we consider the influence 
of charged impurities from the Si02 substrate on the graphene resistivity. This requires 
an estimate of the charged impurity density Uimp in the substrate, which can be obtained 
from a linear fit to the a — n plot at low densities for our graphene samples, prior to the 
addition of the polymerpTTj. This evaluation may be an upper limit since other scatterers 
can also contribute to a linear density dependence of conductivity at low density [15] • How- 
ever, by considering this upper bound, we can at least estimate the maximum contribution 
of charged substrate impurities to graphene resistivity. We obtained an average value of 
nimp ~ 7 X for our samples. Therefore, calculations based on the Boltzmann the- 

ory lead to a maximum contribution of few Ohms for n > The contribution from 

substrate impurities is significantly lower than from electrolyte ions in the polymer matrix. 

The electrolyte ion distribution discussed above is almost temperature-independent since 
the ions are practically frozen below the ice-point of water. Therefore, the phonon contribu- 
tion Pphonon can be extracted from the temperature dependence of the graphene resistivity 
at high-doping. The resistivity versus temperature measurements are shown for Sample 1 
down to 4 K in fig. 3(c). The upturn of resistivity observed at room temperature persists 
down to 4 K, since the resistivity at higher doping remains larger than the resistivity at 
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FIG. 3: (a) Resistivity vs. carrier concentration for Sample 1 (red) and Sample 2 (green), (b) 
Hall mobility vs. carrier concentration for the same 2 samples at T = 295 K. (c) Resistivity vs. 
Temperature at two different densities: n ~ n ~ 2.5 X lO-'^'^/cm^ (blue) (d) 

Carrier concentration vs. applied gate bias; Slope of the linear fit gives an estimate of the gate 
capacitance of the electrolyte gating, C w liiF/cw?. 

lower doping throughout this range of temperature. The temperature-dependent part of 
the resistivity can be fitted by a law at low temperature, up to T ~ 70 K and ~ lOOK 
for n = 2.5 x lO^^/cm^ and 6.2xl0^'^/cm^ respectively. This power-law dependence can be 
associated to a Bloch-Griineisen regime, characterized by a strong suppression of the acous- 
tic phonon scattering rate for T ^ Tbg, where Tbg = 2hvskp/kB is the Bloch-Griineisen 
temperature, Vg the speed of sound in graphene and kp the Fermi momentum [T9]. The den- 
sity of phonons being governed by the Rose- Einstein law, Tbg defines the temperature scale 
below which the acoustic phonon absorption rate vanishes. Besides, the lower the tempera- 
ture, the sharper the Fermi distribution and lower the acoustic phonon emission rate. These 
two factors lead to a complete suppression of the acoustic-phonon induced resistivity in the 
Bloch-Griineisen regime T ^ Tbg- Further, these acoustic phonons are known to be the 
lowest-energy phonons graphene electrons scatter with [TB], which ensures that all phonon 
scattering is suppressed around 4K. Previous measurements down to 20 K and at much lower 
densities n = 2 x 10^^/cm^ to 6xl0^^/cm^ , do not show the Bloch-Griineisen regime [T6] . 
Since Tbg oc ^Jn , the observation of a law to higher temperatures (up to 100 K) in our 
experiment is consistent with theoretical predictions. We also observe a linear regime (or 
non-degenerate regime) between 100 K and 170 K with a slope of PtotiT) ~ 0.13 QK~^, as 
the temperature becomes comparable to Tbg (240 - 420 K), consistent with previous obser- 
vations at lower densities [I6]. Above 200 K, the resistivity becomes a super-linear function 



7 




FIG. 4: (a) Room temperature Hall mobility fitot-ph vs. carrier density for Sample 2. Green 
triangles represent experimental data, shaded blue region represents all possible fitot-ph curves that 
cross the first data point but do not include resistivity contribution from weak scatterers. Orange 
curve is a theoretical fit after including this contribution, (b) Mobility fj^tot-ph vs. carrier density 
for Sample 1 (experimental data in red circles) (c) Resistivity vs. n for Sample 2 (experimental data 
in green traingles). Best fits of ptot = Pion + Pph + Pd + Po to resistivity data: without Fermi velocity 
renormalization (dashed red curve); With electron-electron interaction induced renormalization, for 
/{^lHvK) ~ 0.11 (solid blue curve); By doubling the e-e interaction coupling constant (solid black 
curve). This may reflect the need to include the renormalization from several other interactions as 
discussed in the text. [Inset: Resistivity vs. n for Sample 3] 

of the temperature, indicating that the electrons start to scatter with additional phonons, 
as previously discussed in the literature [HI [161 120] for experiments at lower charge carrier 
densities. Finally, an estimate of Pphonon is obtained as: Pphonon ~ p{295K) — p(4K), which 
is equal to 40 to 47 Q. This almost constant phonon induced resistivity contributes to the 
observed 1/n dependence of mobility. However, even after subtracting Pphonon from the total 
resistivity, the resulting mobility fitot-ph = iX/l^tot — l//^p/i)~^ still shows such dependence 
[see fig. 4(a), 4(b)]. This indicates that other scattering mechanisms are also responsible for 
it, as discussed below. 

We now consider the resistivity induced by defects in the graphene lattice. Strong- 
potential defects such as vacancies and certain adatoms lead to a density-dependent resis- 
tivity of the form, pd = /^,^^rf/[4e^n(ln(i?o^A^))^]) where Rq ~ 1.4 A and are the size and 
density of these defects. Using typical values for nd[T5], this leads to a contribution of the 
order of 1 kQ near the Dirac point. This logarithmic correction leads to a sublinear defect 
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conductivity at high density, and contributes to the decrease in mobihty at high doping. 
However, the increase of nep^ with charge carrier density is too slow to reproduce the mo- 
bihty behavior, as shown in fig. 4. To explain the latter, we thus consider weak scatterers 
inducing a constant resistivity po- This scattering mechanism corresponds to a (5-potential 
of the form, Kca«(^)=Vo5(|r|) [15]. To estimate po, we fitted the theoretical expression 
fitot-ph{n) = l/[ne{pion{n) + Pd{n) + po)] to the corresponding data [fig. 4(a), (b)], giving a 
typical value of ~ 100 Vl. Therefore the mobility analysis shows that the most important 
contributors to the resistivity of our samples at high doping are weak scatterers. 

So far we have discussed the contributions to the graphene resistivity from ions, phonons 
and defects. These contributions are either nearly constant (phonons and weak scatterers) 
or rapidly vanishing with density (charged scatterers and strong-potential defects). It is thus 
surprising to consistently observe a upturn in resistivity in a finite density window near n ~ 
3xl0^^/cm^ [see fig. 3(a), 4(c)]. Note that at higher densities (1.6 xlO^'*/cm^), the resistivity 
decreases, then saturates, as shown for one sample (Sample 3) [see inset of fig. 4(c)]. Below, 
we consider possible corrections to the resistivity terms to model this observed dependence 
on density. At first, we note that phonons and the weak scatterers make up most of the 
graphene resistivity for > 3 x 10^'^/cm^. Therefore, it is natural to consider the upturn 
arising from corrections to these terms. Experimentally, the phonon contribution does not 
increase with density, which makes the weak scatterers the likely cause of the upturn. We 
therefore examine the various factors that determine the resistivity of weak scatterers, in 
order to identify the potential sources responsible for an increased contribution from these 
scatterers. Theoretically, the resistivity contribution of weak scatterers is given by[l5l 119] : 



where Vq is the average impurity potential, uq the density of weak scatterers and vpin) 
is the Fermi velocity. We examine these three factors no, Vq, vp for corrections to resistivity 
induced by weak scatterers. One possible explanation could be the electrochemically-induced 
creation of new defects in graphene, (5no, upon application of gate volatge. However, this 
does not appear very plausible since the onset of upturn is seen in some samples already 
at small gate voltages (e.g. Vg ~ 2V, Ig < 50 pA for Sample 1). At low voltages, the 
resistivity is stable and shows negligible time dependence, precluding the formation of new 
defects which would be a time-dependent process. Furthermore, this resistivity increase is 




(1) 
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reversible with gate voltage and distinguishable from an irreversible increase seen at much 
higher gate voltages {Vg > 10 V), which is likely related to electrochemical processes. A 
different explanation is related to the modification of the local scattering potential Vq at the 
sites where Li^ couples to the carbon lattice to forms complexes. These coupling sites may 
involve a finite density of pre-existing defects on the carbon lattice {e.g. edges, vacancies), 
which could explain the density-dependence of resistivity after the upturn [see fig. 4(c)], as 
these sites get saturated with complex formation. A more interesting source for a resistivity 
upturn is related to the renormalization of Fermi velocity Vp. This renormalization has 
been previously shown by scanning tunneling spectroscopy for graphene on graphite and 
ARPES measurements on epitaxial graphene on silicon carbide substrates in our range of 
densities [5l [21]. A density-dependent renormalization of the Fermi- velocity, vp can lead to 
corrections to an otherwise constant pQ. The Fermi velocity is expected to be renormalized 
by direct electron-electron, Frohlich, electron-phonon and electron-impurity interactions 
\W\ [22H2n]- The direct electron-electron interaction is responsible for an increase of the 
Fermi velocity near the Dirac point, following v — )■ v[l — ln(2a-y/7m)/(4ft:v7i)] but this 
can only partially explain an increase in resistivity in our range of densities [fig. 4(c)]. 
Other factors such as Frohlich and electron-phonon interactions and disorders contribute 
equally to this decrease of Fermi velocity. These factors when considered together can 
potentially explain the magnitude of the observed increase in resistivity [see fig.4(c)] [H [T9l 
[22H2S]- By keeping the concentration of electrolyte ions in the compact layer Cmax around 
10^^/m^ and density of strong-potential defects around 10^^/cm^, it is possible to fit 
the p — n curves of our samples by varying the bare defect resistivity of weak scatterers 
ro around 100 provided that this renormalization of the Fermi velocity is taken into 
account. The upturn is essentially a result of the competition between a decreasing pion and 
an increasing po- Therefore, the charge carrier density corresponding to the onset of the 
upturn increases with the ratio Cmax/fQ- Besides, the decrease in Fermi velocity induced by 
the above-mentioned interactions, slows down at higher densities [n ~ 10^^/cm^), potentially 
leading to a saturation of po [see inset of fig. 4(c)] [22l [211 126H28] . It also follows that for 
samples with a large enough Cmax/fQ, the upturn is expected to be suppressed. Finally, note 
that the Fermi velocity dependence of the phonon-induced resistivity and the renormalization 
of the former do not contradict the fact that the p — T curves [see fig.3(c)] remain almost 
parallel in the linear regime. Due to the screening of the Coulomb interaction between carbon 
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atoms, the sound velocity and the deformation potential decrease with charge carrier density. 
This hmits the influence of a decrease of Fermi velocity on p{T). Before we conclude, note 
that at low densities (n < 5 x the graphene resistivity is already strongly density- 

dependent from Coulomb and strong potential defect scattering. Thus any corrections arising 
from density-dependent renormalization of Fermi velocity in the low density regime are hard 
to discern in transport measurements. 

IV. CONCLUDING REMARKS 

In summary, we have demonstrated high electron densities (6.2xl0^^/cm^) in graphene 
with a polymer electrolyte gate. A Bloch-Griineisen regime was observed between 4 K and 
100 K, a clear sign of large Fermi temperatures. The density-dependence of the mobility 
and resistivity of our samples were analyzed by considering various scattering mechanisms: 
Coulomb scattering from the electrolyte ions, electron-phonon scattering, and electron- 
impurity scattering. Vacancies, cracks and certain adatoms are important scatterers in the 
low density regime. However, weak scatterers are the most important scatterers in our range 
of densities {n > 10^^/cm^), as suggested by the 1/n density dependence of mobility. The 
resistivity versus carrier density graphs, obtained from Hall measurements, show an upturn 
for densities around 3xl0^^/cm^, and possible corrections to resistivity of weak scatterers 
are discussed. While the devices reported in this paper allowed to reach electron densities 
significantly higher than obtained with conventional dielectrics, further improvements are 
needed to explore the physics of graphene in the vicinity of the van Hove singularity. 

V. APPENDIX 

The resistivity induced by the electrolyte ions of the polymer is modeled within the 
framework of Boltzmann theory. Away from the Dirac point, the resistivity is given by p = 
2/[e^Vp{n)D{n)T{n)] where vp is the Fermi velocity, D{n) the density of states, and r the 
scattering relaxation time. We compute the scattering rate induced by a 2D layer of charged 
particles situated at a distance z above the graphene plane by using a screened scattering 
potential, Vc{q,z) = [27re'^e~^^]/[K{q + qTF)]- In a classic Poisson- Boltzmann approach where 
the steric effects are not taken into account, the concentration of ions diverges in the vicinity 
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of the graphene sheet. However, as the gate voltage is apphed, charged ions form a compact 
layer of thickness Ac and concentration Cmax atop graphene[TS]. The total scattering rate 
r^^ is obtained by integrating over the compact layer and the corresponding resistivity reads 
Pionin) = [7ihcmaJion]/[8e\7cnn whcrc lion = /q (1 ' e'^^^'^^") du ^ 0.13, since 

the exponential factor is small for typical values of compact layer thickness (Ac ~ 10 nm) 
and graphene carrier densities {kp ~ 1 nm~^) in our experiment. Note Qs = QTr/kp in the 
integral. 
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